Mean field approach to flavor susceptibilities with a vector interaction 
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We show that flavor diagonal and off-diagonal susceptibilities of light quarks at vanishing chemical 
potential can be calculated consistently assuming the baryon density and isospin density dependence 
of QCD to be expressed by a vector-isoscalar and a vector-isovector coupling, respectively. At the 
mean field level, their expression depends only on the effective medium-dependent couplings and 
quark thermodynamic potential. The strength of the couplings can be then estimated from the 
model using susceptibilities calculated in lattice QCD as an input. 
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, . I. INTRODUCTION 

■ The recent observation of the large elliptic flow at RHIC P, 0| together with the high jet energy loss, led to the 
conjecture that the matter at RHIC is strongly coupled, a nearly perfect fluid. On the other hand, Monte Carlo 
simulations of QCD (Quantum Chromodynamics) on a discrete lattice (lattice QCD) I3| suggest a quasiparticle 
picture of the Quark Gluon Plasma (QGP), at least for the quarks. In fact, both flavor-off-diagonal susceptibilities [3| 

■ and higher order baryon number susceptibilities Q are consistent with vanishing correlations for teinperatures right 

■ above the transition, T > l.2Tc [1], and appear to be well described by quasiparticle approaches In addition, 
I ] effective quark models with a Polyakov loop potential (see, for example, refs. d, 0) proved to be quite successful in 

^ reproducing the thermodynamics of strongly interacting matter, and an even better agreement was recently found 
^ with a (2+l)-flavor Polyakov-quark-meson model JT|j in the mean-field approximation. Based on these observations, 

^ ^ it was proposed in ref. [(| that the large elliptic flow at RHIC could be described by single-particle dynamics with 
a repulsive interaction. The same interaction would also explain the (~ 15%) deviation from the Stefan-Boltzmann 

psj limit of the pressure calculated in lattice QCD. 

J> A repulsive mean-field of the vector type was also invoked as a possible explanation for the discordant results of 

00 lattice QCD in refs. that appear to disfavor the existence of a critical point in the QCD phase diagram. The 
authors of refs. [ll| observed that the region of quark masses where the transition is of the first order (for quark 
masses smaller than the physical ones), tends to shrink for small positive values of the chemical potential ^. This is 
in contrast with model studies that support the existence of a critical point. In these, the first order region has to 
' expand with increasing fi, so that the physical quark mass point coincides with the critical line at some finite value of 
T and /i. As was proposed in ref. [l^ , a sufficiently strong vector coupling may account for the shrinkage of the first 
order region, eventually resultin g in two critical points for a given range of (small) quark masses and in one critical 
point for physical quark masses [l^ ■ 

Inspired by these ideas, we propose a scheme to estimate the strength of the vector interaction from lattice QCD 
data on flavor diagonal and off-diagonal susceptibilities. We will confine ourselves to the light quark doublet and 
we will assume that, at high temperature, the dependence on the net-baryon density and the isospin density can 
Ci ■ be described by mean field vector-isoscalar and vector-isovector interactions, respectively. This is somehow in the 
same spirit of ref. where it was shown that lattice data on baryon number susceptibility may be described in 
terms of a temperature-dependent repulsive vector coupling in the Nambu Jona-Lasinio (NJL) model. The vector- 
isovector coupling, needs also to be taken into account as its relative strength, with respect to the vector-isoscalar 
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coupling, has a strong impact on the flavor off-diagonal susceptibilities [15[- In our analysis, we do not make any 
specific modeling of the remnant part of the interaction (that is certainly present in QCD). Its effect is taken into 
account by allowing the vector couplings, as well as the quark thermodynamic potential^, to be medium-dependent 
parameters. Within these assumptions, we will consistently obtain a very simple expression for the flavor diagonal and 
off-diagonal susceptibilities at zero chemical potential. The result can be then inverted and the difference between the 
vector-isoscalar and the vector-isovector coupling can be extracted directly from lattice data using the susceptibilities 
evaluated in [l6j as an input. The strength of both couplings can be then estimated considering two extreme cases that 



As it will be clarified in the following, here we use the expression "quark thermodynamic potential" to indicate the thermodynamic 
potential of the quark sector at vanishing vector coupling. 
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will allow us to obtain lower and upper bounds for these quantities as a function of the temperature. As we will show, 
the vector-isoscalar coupling is very large (repulsive interaction) immediately below and reaches quickly much lower 
values above Tc, where it basically coincides with the vector-isovector coupling. Below Tc also the vector-isovector 
coupling exhibits a similar behavior but remains always smaller than the vector-isoscalar coupling. 

This note is organized as follows: In Sec.|TT]we briefly review the basic mean field results for a massive vector field 
interaction and we introduce the partition function of the model. In Sec. Illll we calculate net quark number densities 
and susceptibilities and in Sec. II VI we extract the value of the vector couplings from lattice QCD data. 



II. THE MODEL 



In this section, we will introduce the effective partition function of our model. To start, we will recall a few basic 
results of mean-field theory in presence of a vector-isoscalar coupling gy and a vector-isovector coupling gy. The 
Lagrangian of a system of quarks with such massive vector fields is of the form 

£ - * [i^d^ - m) VI/ - -^F^,F^^- - glr'^^-^V^ + y ^^V^ - gl^^T^B^ + ^B^B^ , (2.1) 

where vj/ = (u, d) is the quark field, ra is the quark bare mass, r = (r^, r^, r'^) are the Pauli matrices and 77^, "Hv are 
the bare masses of the vector-isoscalar field and the vector-isovector field B^, respectively. In the Euler-Lagrange 
formalism we can calculate the equation of motion for the static vector fields. By taking the average we obtain for 
the 0-th component of V"^ = (Vb, V) and = {Bq, B) 

m = ; Pq ^ {^^^) (2.2) 

The average of the spatial parts (V) and (B) vanishes because of the rotational symmetry. The same is true for {Bq) 
and (-Bq), because of isospin symmetry. In mean field approximation, the eigenvalues of the energy of the quarks turn 
out to be (see for example jla, llTj) 

^^^(P) = ^p^ + ^pi±Vp^+m^ (2.3) 

vi vi 
Edip) = ^P^-^Pi±V^^T^ , 

with the energy for anti-u(d) quarks being — _E~|.^^(p), according to the ordering prescription in Following the 
standard procedure, one can write down the partition function for this simple system. For our discussion, we will 
confine ourselves to the light flavor doublet (u and d quarks) and we will adopt classical (Boltzmann) statistics. After 
normal ordering one obtains 



^ ^cosh(^ -" -^-- j$„ + cosh(^ -" -v-^^ ■ + + (2.4) 

where fj.^ and pd are the chemical potentials for the u and d quarks and — and —^dT are the corresponding 
thermodynamic potentials for Pu = Pd = 0. In this simple case and $d are just the logarithm of the partition 
function of a gas of free quarks with mass equal to the bare mass m. At the mean field level, the couplings Af^ = 
9v'^/Vs and Ay = gv'^/ilv are constant. In Eq. (12.41) . the terms \yp^/2T and Xypf/2T are sometimes referred to as 
rearrangement terms. Additional terms of this kind appear naturally to account for the energy balance and restore 
thermodynamical consistency in mean field approximation [18] . Our analysis will be performed in the infinite volume 
limit. From now on, to simplify the notation, we will omit the system volume V in the following equations. 

The Eq. (|2.4p is indeed an oversimplified partition function to describe a complicated theory such as QCD. Even if 
we want to confine ourselves to a purely phenomenological description through effective models, we need, at least, to 
include other kinds of fields and interactions. For some purposes, it may be sufficient to add a scalar field (as in the 
Walecka model [l3|) or a scalar and a pseudoscalar meson (as in the linear sigma model with quarks [l^) with their 
interaction terms. For our discussion, however, we will not need to know the exact nature of these interactions. We 
rather assume that the vector-like part of the interaction in the QCD Lagrangian can be isolated and treated as a 
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mean field. In the quark sector, the effect of the rest of the interaction is then assumed to be relegated to the vector 
couplings and the quark thermodynamic potential. Starting from Eq. (|2.4[) we will consider a partition function of 
the form 

= cosh f ^^~^v*P^.-K*PA ^ ,,3, f ,.-K*P^. + K*PA + 1 + l^pf + R . (2.5) 



In Eq. (|2.5p . the couplings Xy* and Ay* are now medium-dependent quantities that may embody radiative corrections 
arising from additional pieces in the Lagrangian in Eq. (|2.ip . Similarly, also and can have a non-trivial 
dependence on the medium. Any kind of additional interaction, besides affecting directly the couplings and the quark 
thermodynamic potential, will in general also result in additional pieces in the partition function. For instance, in 
presence of (pseudo)scalar excitations (besides the appearance of a thermal quark mass implicitly embodied in $* and 
^J), the Eq. (|2.5p should contain terms such as the partition function of the (pseudo)scalar quasiparticles plus any 
rearrangement term needed for the thermodynamic consistency of the modeP. These contributions are represented 
by the term R in Eq. (|2.5p . As we will show below, this term does not enter in the expression for the densities and 
the susceptibilities, but must be taken into account to evaluate the pressure. 

Equation ()2.5p exhibits the whole explicit dependence on pq, pi, ^u, Md, Ay* and Ay* of the partition function. 
In general, the quark thermodynamic potentials, the couplings Ay*, Ay* and R will depend on T and on a set 
{Hi (/Xu, /Xd, r)} of additional quantities stemming from the remnant part of the interaction. In a thermal environment, 
these terms could be the scalar condensate or, as in ref. [13, [131 the average fluctuations of the scalar mesonic fields 
and so on. Formally, our partition function can be written as 

In Z ^j.d,T) = In Z [^u, Md, T, pq (/iu, Md, T") , Pi (Mu, Md, T") , {n* (^u, Md, T")}] ■ (2.6) 

Following the general prescription in ref. ^211] we then impose the thermodynamical consistency through the condition 
for the stationarity 



d\nZ 



dpq 



dl.Z 



dpi 



and 9lnZ 



= Vi . (2.7) 



Eqs. (|2.7I) need to be fulfilled by any model to recover the standard connection between statistical mechanics and 
thermodynamics^ (for a detailed discussion see ref. [2l|). In other perturbative approaches, such as the Optimized 
Perturbation Theory an analogous condition is imposed by the principle of minimal sensitivity [1^ to determine the 
value of the thermal mass (see for example refs. [24| and references therein). In the following we will make use of 
these two equations without making any further assumption on the functional form of Ay*, Ay*, $*, $J and R. To 
simplify the notation we recast the Eq. (|2.5p in a more convenient way by using the isospin symmetry 



$ : = $* ^ $ (2.8) 



and the shorthand notations 



leading to 



\S * \V * II II 

= — ; >^v = — ; ttu = — ; ttd = Y ' (2.9) 



InZ ^ cosh (au - XsPq - Kpi) $ + cosh (ad - XsPq + Kpi) * + ^^sPq + ^K'Pi + ^ ■ (2-10) 



By further introducing the chemical potential /iq for the net quark-number A^q = + Nd and the chemical potential 
pi for the isospin'' I = 2I3 — — Ad one defines 

"q = Y " ^ + ' "I ^ ^ = ^ ("u - ad) , (2.11) 



^ These rearrangement terms are always present in density-dependent field theories where they appear in the self-energy corrections [2Cl| . 
In some cases they may be interpreted as a consequence of the medium-dependent energy of the system in absence of quasiparticle 
excitations [2l| . 

^ The first equation in Eqs. 1 12.711 should not be confused with the differential of the Gibbs free-energy G at constant temperature and 
pressure: G = fiN, (where A'^ is the number of particles and fi the chemical potential). Here pq is not a natural variable of the partition 
function (thermodynamic potential), and the stationarity condition has nothing to do with a variation of the physical quantities in the 
system. The same arguments hold for pi. 

* Note that we count isospin in units of 1, instead of 1/2. 
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and the partition function finally reads 

In Z = 2 cosh (aq - A^pq) cosh (ai - A„pi) $ + ^ A^Pq + ^Kpi + R . (2.12) 



III. DENSITIES AND SUSCEPTIBILITIES 



The net quark-number density pq and the isospin density pi can be derived from the partition function in Eq. (j2.12p 
by performing the partial derivative with respect to aq and ai, respectively. The dependence on aq, ai of pq and pi, 
and the implicit dependence on aq, ai of As, $ and R can be ignored due to Eqs. (|2.7p . This results in the implicit 
expressions 

9 In Z 

= Pq = 2sinh(aq - A>,Pq)cosh(ai - At,pi)$ , (3.1) 

= Pi = 2 cosh (aq - A^pq) sinh (aj - A^,pi) $ . 

We note that, thanks to the rearrangement terms jAsPq and ^Xypf, the first two equations in Eqs. (|2.7|) are naturally 
fulfilled by the solutions in Eq. p.ip 



dlnZ 



dpq 
dlnZ 



dpi 



= -2As sinh (aq - A^pq) cosh (aj - At,pi) $ + A>,Pq = , (3.2) 
~ —2Xy cosh (aq — AsPq) sinh (aj — At,pi) $ + X^pi = . 



Similarly, the other rearrangement term(s) in R must guarantee the validity of the third relation in Eq. (|2.7p . From 
Eq. (|2.10p . one can also find the u and d quark net number-densities 

Pu = sinh (au - A^Pq - XvPi) $ ; Pd = sinh (ad - A^pq -I- X^pi) $ . (3.3) 

Because of charge conjugation- and isospin-symmetry all densities defined in Eqs. (j3.ip and p.3p vanish for aq = aj = 
0. More generally, any odd derivative of the partition function vanishes. 

Now that we have the densities, we can go ahead and calculate the susceptibilities. These are obtained by further 
differentiating with respect to the chemical potential. Introducing the notation a = (aq, ai) we define: 

1 dpa 2 > dpa dXs 



Xq{a,T) = -77-^ = -cosh(aq - AsPq)cosh(ai - A„pi)$ 1 - As—-^ - Pq 

1 OUq 1 \ OUq OUq 

2 / dp d\ \ 
- — sinh(aq - AsPq)sinh(ai - A„pi)<I> {^Vq^ '^^'da') 

+ 2 sinh (aq - AsPq) cosh (ai - Ai,pi) - — ; 

oaq 

1 dp 2 / dp d\ \ 
XI (a, T) = — ^ = -—sinh (aq - AsPq) sinh (ai - Aupi)$ fAs^ +Pq^j 

2 / dp d\ \ 
+ ^cosh(aq - AsPq) cosh(ai - A„pi) $ f 1 - -^u^ - '^^'dc^/ 

2 (9$ 
+ — cosh (aq - AsPq) sinh (ai - A-upi) — . (3.4) 

For a = these two equations reduce to: 

2$ 



Xq = 



T(H-2As$) 



2$ 

= T(1-|-2A,$) ' ^^-^^ 

where we have used the shorthand Xq = Xq(0,r) and the same for xi- The terms involving the derivatives of As, A^, 
and $ with respect to aq and ai in Eqs. p.4p do not contribute to the susceptibilities for a = 0. Unfortunately, this 
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is not the case for higher order susceptibihties. Aheady at the 4-th order there are non vanishing terms involving the 
second derivative^ of Ag and A^,. To evaluate these terms, one then needs to introduce further assumptions in the 
model, but this goes beyond the scope of our discussion. 

The diagonal and off-diagonal flavor susceptibilities Xuu and Xud are given by Xuuicx-,T) = {l/T)dpu/dau and 
Xud {(^,T) — {l/T)dpu/dad which can be conveniently obtained using (see Eqs. (|2.1ip ') 



d_ _ 1 f_d_ d \ 



dau 2 \daq dai J 

_d_ _ 1 f_d d_\ 

dad 2 \dan dai J 



(3.6) 



Therefore, 



1 92 In z 1 / 92 ^2 



For a = the mixed derivatives on the first equation vanishes, and therefore: 

Xq + Xl . 



X 



uu 
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Xud = ^^^^ . (3.8) 

Note that, for a = 0, we have Xuu = Xdd because of isospin symmetry. From Eq. p.Sp and Eq. (13.81) one has: 

$[1 + (A, + A,)$] 



X 



UU 



r(l + 2A,$)(l + 2A„$) 



_ (A. - A,) 

Xud - r(i + 2A,$)(l-f 2A„$) • ^ ' 

It is worth to remark that this result, together with Eqs. p.Sp . does not formally depend, within our model assump- 
tions, on the possible presence of other kinds of effective interaction. Their effect is implicitly taken into account by 
the non-trivial dependence on the medium of the functions $, A^ and At,. Particularly, the presence of other (scalar) 
mesons would affect directly the value of $ (and therefore of the susceptibilities), for instance through a temperature 
dependent effective quark mass. There are, however, other mesonic effects that are not included in our scheme. As 
shown by Eq. p.9p . the off-diagonal susceptibility Xud vanishes if the isovector and isoscalar couplings coincide (or 
vanish) . This cannot account for the contribution of pions described in ref. i25t] , which exist also in the absence of a 
vector interaction. A more detailed discussion may be found e.g. in refs. [^H^- Inverting the Eqs. p.9p one finds: 

In principle, the vector couplings could be then estimated from Eqs. (|3.10p using lattice QCD results on flavor 
susceptibilities as an input. It is, in fact, tempting to identify the pressure for a = with 2$r. However, one should 
bear in mind that we are actually missing the contribution RT coming from the unknown term R in Eq. (j2.10 ). As 
an example, in the presence of additional scalar and pseudoscalar mesons as in the linear sigma model of ref. [22 1, the 
factor R will account for the pressure generated by the mesons, plus other terms that must be introduced in order 
to avoid double counting of the fundamental degrees of freedom. Unlike Xq and xij the information carried by is 
therefore incomplete and cannot be immediately related to the pressure of the quarks. In the next section, however, 
we will show that is possible to explore two limiting cases and provide upper and lower bounds for A^ and A„. 



For symmetry reasons, one expects <&, and A„ to depend on even powers of Oq and aj, resulting in vanishing first derivatives for 
q: = 0. 
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IV. VECTOR COUPLINGS FROM LATTICE QCD 



In this section, we will provide an estimate for the temperature-dependent vector cou plin gs As and starting 
from Eqs. (|3.10p . To do that we will use, as an input, the susceptibilities evaluated in ref. [ig in a two- flavor lattice 
simulation with bare quark mass m/T = 0.4 at a = 0. In ref. [16] it was shown that Xq and Xi (and then, due to 
Eq. (|3.8I) . also Xuu) rise steeply around Tc and then reach a plateau at approximately 80% of the Stefan Boltzmann 
(SB) limit. The off-diagonal susceptibility Xudj instead, starts from negative values {xi is a bit larger than Xq at 
T <Tc) and then approaches zero at T > T^. The vanishing correlations between u and d quarks do not imply that 
there are no interactions. The deviation from the SB limit indicates that they are actually still important at T ^ 2Tc- 
From the second equation in Eqs. (13. 9|) one sees that Xud = if At, = Ag. The deviation from the SB limit can be 
then due to vector-type interactions, and/or to non- vector interactions reflected directly in the quark thermodynamic 
potential $ (see Eqs. p.5p ). 

From Eqs. p.lOp . we note that the difference (As — A„) does not depend on $ 

As-A. = lf--1) , (4.1) 

T \Xq XI J 

and can be evaluated using lattice data. The result is shown in Fig. [T]i. As one can see, the dimensionless quantity 




FIG. 1: Left panel: the difference (As — \v)T^ extracted from lattice QCD data [l6l ] on flavor diagonal and off'-diagonal 
susceptibilities using the formula in Eq. (|4.ip . Right panel: The assumed upper bound (<I>*Vr^) and lower bound (<E>™/r^) 
for the function $/r^. 



(As — \v)T^ is large for T/Tc < 1 and vanishes very quickly above the critical temperature. The vector-isoscalar 
coupling As is therefore larger or equal (for T > Tc) to A^. This was somehow expected as in the chirally broken 
phase the coupling of the vector-isoscalar field (the uj meson) with the nucleons is empirically three times larger than 
the corresponding vector-isovector (the p meson) coupling, whereas in the symmetric phase they are expected to be 
equal (see, for example, the discussion in ref. [lq|). 

As already mentioned, a unique estimate for the couplings cannot be obtained because the function $ in Eqs. p.lOp 
cannot be extracted from lattice data without introducing additional assumptions to quantify the effect of the unknown 
contribution to the pressure RT. From Eq. (|4.ip it is clear, however, that the lower bound for As, A™ corresponds to 
vanishing At, . In this case A™ is simply given by the Eq. (|4.ip with At, = and the second relation in Eq. (I3.10p gives 



(4.2) 
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where the superscript m means that this quantity corresponds to the lower bound A^. By plugging xi from lattice 
QCD into the Eq. (I4.2p one can now evaluate the function The result (divided by T^) is the continuous curve 
with open circles in Fig. [TJd. In this scenario the coupling goes to zero very rapidly above Tc (since A„ = the 
curve in Fig. [T^ corresponds to X™T^) and only interactions which are not of the vector type, such as (pseudo)scalar 
interactions, are left. The latter entirely account for the aforementioned 80% deviation from the Stefan Boltzmann 
limit of the susceptibilities. As a result, the thermodynamic potential is much smaller than the corresponding SB 
limit for massless particles 12T^/7t^. 

The opposite extreme case, instead, is realized by imposing that the only interaction responsible for such a deviation 
is of the vector type. As one can see from the first equation in Eqs. (|3.10p . A^ is maximal when $ is maximal. To 
estimate the upper bound of As and A„ we therefore replace $ — )• in Eqs. (|3.10p . where = and A is 

a constant value chosen is such a way that ^^'^ (2Tc) = SB limit (see Fig. [TId). This will give us the upper bounds 
A^^ and A^"^ corresponding to the case where the whole deviation from the SB limit at T > 2Tc is due to the vector 
interactions^. This is quite an overestimate, as quarks are actually expected to have a thermal mass as found e.g. 
in the hard thermal loop expansion (see, for example, refs. [IBl)- By construction, therefore, our estimates leave the 




U I I I I I I 

0.8 1 1.2 1.4 1.6 1.8 2 

T/T 



FIG. 2: The vector-isoscalar coupling As and the vector-isovector coupling Xy (multiplied by T^) as a function of T/Tc. Shaded 
areas are drawn between the upper and the lower bounds representing our estimated uncertainty. 



possibility of an important vector coupling at high temperature (as was proposed in ref. Q ) open. To obtain tighter 
constraints, other (independent) flavor-related observables should be considered. One possibility, could be higher 
order susceptibilities. This, however, requires the introduction of additional assumptions as mentioned in Sec. IIIII 
The admissible values for As and A^ are represented by the shaded areas in Fig. [21 As one can see, the dimensionless 
quantity AsT^ is large (repulsive interaction) for T/Tc < 1 and reaches quickly much lower values above the critical 
temperature. A similar behavior is shown also by Xy that, however, is smaller than As at low T. 

As mentioned in the introduction, the effect of a vector-isoscalar coupling was also accounted for in a recent work [12 1 
to describe the shrinkage of the first order region in the mudj^s plane (where rriud is the light quark mass and mg 



Note that if there are only interactions of the vector type the term R in Eq. 1 12. 121 1 is zero and 2<3?T is the actual pressure of the quarks. 
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the strange quark mass) at small chemical potential /x. This effect was observed in lattice QCD [llj and led to the 
conclusion that the existence of a critical point (at least at small /i/T) is unlikely. The author of ref. [l^l was able to 
reproduce this behavior in the Polyakov loop extended NJL model (PNJL) by introducing a repulsive vector- vector 
interaction of the isoscalar type. For a sufficiently strong vector coupling'' (Ay ^ 0.8A, where A is the scalar coupling) 
the first order region initially tends to shrink as fj, is increased to small positive values, and, after reaching its minimum 
size, it tends to expand as /i is further increased. 

The method proposed in this note could be, in principle, used to extract the value of the vector-isoscalar coupling 
in coincidence with the critical line T = Tc and m'^^ = m^, where and are the chiral critical masses. To do 
that, however, further ingredients are needed. First, our treatment should be extended to three flavors. Second, also 
the scalar coupling should be extracted from lattice data in a consistent fashion. In effective quarks models such as 
the NJL model, in fact, the scalar coupling A depends on the cutoff and the value commonly used in literature [27l.[28j 
ranges from A 5 GeV~^ to A ~ 10 GeV~^. To provide a rough comparison (with all the aforementioned caveats in 
mind) from our analysis we obtain, assuming Tc = 0.17 GeV, 4 GeV^^ ^ ^ 19 GeV~^. At this level, therefore, 
no conclusions could be drawn in this respect. 

V. SUMMARY 

We proposed a method to estimate the strength of the vector interaction from lattice QCD data on flavor diagonal 
and off-diagonal susceptibilities. Our fundamental assumption is that the baryon density and isospin density depen- 
dence of QCD can be described as a mean field by a vector-isoscalar and a vector-isovector coupling, respectively. 
Imposing thermodynamical consistency conditions, we showed that in our framework flavor susceptibilities can be 
parametrized by the medium-dependent effective couplings Ag (vector-isoscalar), A^, (vector-isovector) and by the 
quark thermodynamic potential <&. The effect of other kinds of interaction, like e.g. scalar or pseudoscalar, affects the 
value of these quantities, but not the form of the result. Inverting our relations, we have shown that the difference 
(As — Xy) depends only on Xq and xi (or, equivalently, Xuu and Xud) and can be then extracted from the model using 
the susceptibilities from lattice QCD |Tq| as an input. We found that (As — \y) is quite large below Tc and approaches 
zero very rapidly above Tc- 

Unlike their difference, the value of the vector couplings alone depends on the details of the remnant part of the 
interaction. Therefore, they cannot be determined directly from lattice QCD without introducing further assumptions. 
We have then estimated upper and lower bounds for As and A„ (using again lattice data as an input) considering 
two quite extreme scenarios. The first (yielding the lower bounds), where the effect of the non- vector part of the 
interaction (such as scalar or pseudoscalar) is maximum and the second (yielding the upper bounds) assuming that 
the only interaction left at T > 2Tc is of the vector type. Immediately below Tc, our results suggest a strong As 
that rapidly approaches lower values at T > Tc- A similar behavior is observed for Xy that, however, remains always 
smaller than As at low temperature. 

In future works this method should be refined in order to provide tighter bounds for the possible values of the 
couplings. This could be done considering other flavor-dependent observables from QCD (such as higher order 
susceptibilities) and/or motivating realistic assumptions to estimate the bounds. Our analysis should be also extended 
to three flavors. This could be helpful to understand if the strength of the vector coupling As is large enough to produce 
a change of sign of the curvature of the critical surface in the mud, "^g, fi space as proposed in ref. [12|] . 
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Notice that the author of ref. [T^ uses a different definition for the vector coupling Gy = 
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